AD-A245  876 

.'Hill  ....  _  ^  *  W 


NAVAL  FUM  GRADUATE  SCHOOL 

Monterey,  California 


DTIC 


elects; 

FEB141S92|’ 


0 


THESIS 


HIGH-GAIN,  HIGH-POWER  FREE  ELECTRON  LASERS 

by 

Jung-Hyun  Park 

June,  1991 

Thesis  Advisor. 

W.  B.  Colson 

Approved  for  public  release;  distribution  is  unlimited 


92-03537 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

1b  RESTRiaiVE  MARKINGS 

2a  SECURITY  CLASSIFICATION  AUTHORITY 

3  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution  is  unlimited. 

2b  DECLASSIFICATION/DOWNGRADING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6a  NAME  OF  PERFORMING  ORGANIZATION 
Naval  Postgraduate  School 

6b  OFFICE  SYMBOL 
(If  applicable) 

33 

7a  NAME  OF  MONITORING  ORGANIZATION  1 

Naval  Postgraduate  School 

6c  ADDRESS  (C/ty,  State^  and  ZIP  Code) 

Monterey,  CA  93943-5000 

i 

7b  ADDRESS(C/fy,  State.  andZIPCode) 

Monterey,  CA  93943-SOOO 

8a  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

8b  OFFICE  SYMBOL 
(If  applicable)  j 

9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

8c  ADDRESS  (Crty,  State,  and  ZIP  Code) 

10  SOURCE  OF  FUNDING  NUMBERS 

Program  tiemeni  No 

Projea  No 

1  dSk  Nc 

Work  Unit  Acces^on 

Number 

11  fULt  (Include  Security  Classification) 

HIGH  GAIN,  HIGH  POWER  FREE  ELECTRON  LASERS 

12  PERSONAL  AUTHOR(S)  PARK,  JUNG  HYUN 


13a  TYPE  OF  REPORT 

13b  TIME  COVERED 

14  DATE  OF  REPORT  (year,  month,  day) 

15  PAGE  COUNT 

Master’s  Thesis 

From  To 

1991  JUNE 

42 

16  SUPPLEMENTARY  NOTATION 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reOect  the  oflicial  policy  or  position  of  the  Department  of  Defense  or  the  U.S. 
Government _ _  _  _ ' 


17  COSATI  CODES 

1 8  SUBJECT  TERMS  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

FIELD 

GROUP 

SUBGROUP 

FEL,  undulator  .gain  spectrum,  gain  evolution 

19  ABSTRACT  (continue  on  reverse  if  neceuary  and  identify  by  blotk  number) 


The  LLNL  Paladin  FEL  experiment  is  shown  to  exhibit  clear  and  dramatic  effects  governed  by 
the  electron  beam  velocity  distribution  for  the  first  time.  The  FEL  integrsd  equation  is  used  to  show 
that  there  is  significant  broadening  of  the  gain  spectrum  due  to  the  Gaussian  velocity  distribution, 
and  also  shows  a  plateau  in  the  gain  evolution  along  the  undulator  due  to  a  triangular-shaped 
velocity  distribution. 

The  gain  spectra  and  power  evolution  from  simple,  single-mode  simulations  are  compared  to 
the  LLNL  ELF  experiments.  The  microwave  power  evolution  along  the  undulator  is  compared  as 
well  for  both  the  tapered  and  untapered  undulators.  In  all  cases,  the  agreement  is  found  to  be  good. 


20  DlSTRlBUTtON'AVAILABlLlTY  OF  ABSTRACT 

n  UNCtASSifitDilNl'MlItU  Q 

21  ABSTRACT  SECURITY  CLASSIFICATION 

unclabaifieti 

22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

W  H  ('olsoii  j 

22b  TELEPHONE  (Include  Area  code) 

(4081 646  276.5 

DD  FORM  1473,  84  MAR  83  APK edition  may  be  u^ed  until  exhausted  ClASSiMCATIQN  Qt  THIS  PACt 

All  other  edilinna  are  obsolete  U  NCl.ASSI  FI  K  f) 


I 


Approved  for  public  release:  distribution  is  unlimited. 


HIGH-GAIN,  HIGH-POWER  FREE  ELECTRON  LASERS 

by 

Jung-Hyun  Park 

Captain,  Republic  of  Korea  Army 
B.  S.,  Korea  Mibtary  Academy,  1985 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  PHYSICS 
from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1991 


Author: 


Approved  by: 


H. 


Jung-Hyun  Park 

U)  ■  'S-  e-^ 


William  B.  Colson,  Thesis  Advisor 


Department  of  Physics 


ii 


ABSTRACT 


The  LLNL  Paladin  PEL  experiment  is  shown  to  exhibit  clear  and  dramatic 
effects  governed  by  the  electron  beam  velocity  distribution  for  the  first  time.  The 
PEL  integral  equation  is  used  to  show  that  there  is  significant  broadening  of  the 
gain  spectrum  due  to  the  Gaussian  velocity  distribution,  and  also  shows  a  plateau 
in  the  gain  evolution  along  the  undulator  due  to  a  triangular-shaped  velocity 
distribution. 

The  gain  spectra  and  power  evolution  from  simple,  single-mode  simulations 
are  compared  to  the  LLNL  ELP  experiments.  The  microwave  power  evolution 
along  the  undulator  is  compared  as  well  for  both  the  tapered  and  untapered 
undulators.  In  all  cases,  the  agreement  is  found  to  be  good. 
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I.  INTRODUCTION 


In  a  Free  Electron  Laser  (FEL),  a  relativistic  high-quality  beam  of 
electrons,  passing  through  an  undulator  magnet,  generates  a  powerful  and 
coherent  electromagnetic  wave  [1].  Coherence  may  develop  after  many 
passages  in  the  undulator,  when  the  device  operates  as  an  oscillator,  or  in 
one  passage  in  the  amplifier  configuration.  The  choice  between  an  oscillator 
configuration  and  a  amplifier  configuration  is  generally  determined  by  the 
accelerator  that  drives  the  FEL;  specifically  by  the  pulse  format  and  peak 
current  of  the  accelerator.  Induction  linear  accelerators  produce  short 
pulses,  usually  10-100  ns  long,  that  are  widely  spread  in  time  with  high 
peak  current.  The  pulses  are  too  far  apart  to  be  synchronized  with  optical 
pulses  in  a  resonator.  To  be  interesting,  the  amplifier  must  have  high  gain 
in  a  single  pass  of  the  light  through  the  amplifier.  The  high  peak  current  of 
an  induction  linac  permits  both  high  gain  and  high  efficiency,  the  fraction  of 
electron  beam  energy  converted  to  light. 

The  evolution  of  coherent  electron  bunching  is  the  key  element  to  any 
FEL  gain  mechanism.  In  a  high-gain  amplifier,  like  the  Paladin  FEL  [2-4] 
at  Lawrence  Livermore  National  Laboratory  (LLNL),  maintaining  coherence 
of  the  electron  bunch  over  a  significant  interaction  length  imposes  important 
restrictions  on  the  initial  electron  beam  quality.  The  design  of  an  high  gain 
FEL  amplifier  requires  high  beam  current  and  a  long  undulator;  both 
features  tend  to  allow  gain  to  be  easily  degraded  by  poor  beam  quality. 


1 


Paladin  is  the  first  short-wavelength  PEL  to  operate  in  the  high  current 
regime  where  the  growth  mechanism  is  dominated  by  the  longitudinal 
velocity  distribution  of  the  electron  beam.  It  is  observed  that  the  growth  of 
the  optical  power  along  the  15  m  Paladin  undulator  is  not  exactly 
exponential,  but  is  determined  by  the  shape  of  electron  distribution 
function.  Similarly,  the  gain  spectrum  measured  in  the  earlier  5  m  long 
Paladin  amplifier  shows  a  dramatically  broader  and  lower  shape  due  to  the 
electron  beam  distribution.  The  PEL  integral  equation  [5-7]  is  used  to 
calculate  the  optical  field  evolution  in  the  presence  of  an  arbitrary  electron 
velocity  distribution. 

The  LLNL  Electron  Laser  Pacility  (ELP)  PEL  is  a  single-pass,  high- 
gain,  high-power  amplifier  operating  in  the  microwave  regime  at  34.6  GHz 
frequency  [8-11].  The  characteristics  of  this  PEL  have  been  extensively 
studied  with  large  simulations  including  many  transverse  modes, 
diffraction,  beam  emittance,  and  betatron  focusing.  A  relatively  simple, 
single-mode  model  of  the  PEL  interaction  using  only  the  electron  pendulum 
and  optical  wave  equations  can  describe  the  most  important  features  of  the 
ELF  PEL.  The  waveguide,  used  in  the  experiment  to  confine  the  natural 
diffractive  spreading  of  the  microwave  beam,  can  be  simply  modeled  by  a 
filling  factor  and  a  shift  in  resonance. 

The  most  significant  contributions  to  the  analysis  of  a  high-gain,  high- 
power  PEL  are  summarized  below. 

•  Gain  spectrum  broadening  of  the  5  m  Paladin  experiment  due  to 
poor  beam  quality  are  examined  by  PEL  integral  equation. 
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Plateau  in  the  15  m  Paladin  PEL  gain  evolution  is  explored  by 
the  velocity  distribution  functions. 

The  asymmetrical  shape  of  high-gain  ELF  PEL  is  examined  by 
simple  PEL  equations. 

Gain  evolutions  with  both  untapered  and  tapered  ELF  PEL  are 
explored  using  simple  PEL  equations. 


II.  BASIC  FEL  THEORY 


A  ELECTRONS  IN  THE  UNDULATOR 

In  order  to  discuss  the  principles  of  the  FEL,  it  is  important  to  consider 
the  individual  electron  trajectories  in  the  undulator  field  alone,  without  the 
presence  of  light.  The  interaction  between  the  negatively  charged  electrons 
is  negligible  in  tj^pical  FELs  [1],  and  the  beam  is  focused  externally  before 
entering  the  undulator.  Under  the  influence  of  the  undulator  magnetic 
field,  the  electrons  will  deflect  periodically  in  the  transverse  direction,  or 
"wiggle"  as  they  travel  along  the  optical  mode  axis.  The  ideal  form  of  the 
on-axis  helical  undulator  field  atx  =y  =  0  is 

=  B[cos(^o2’ ,  for  0<2<L  ,  (1) 

where  B  is  the  peak  undulator  field  strength,  kQ  =  2tiI'Kq  is  the  undulator 
wavenumber  associated  with  the  wavelength  z  is  the  distance  along  the 
longitudinal  axis  of  the  undulator,  L  =  is  the  undulator  length,  and  N 
is  the  number  of  undulator  periods.  The  Lorentz  force  equations  [12] 
become 

,  and  -^=0  ,  (2) 

dt  me  dt 

where  is  the  Lorentz  factor,  m  is  the  electron  rest  mass, 

e  -  \e\  is  the  electron  charge  magnitude,  c  is  the  speed  of  light, 
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P  =  t!*/C  =  VK  .  is  the  perpendicular  component  which  refers  to  the 
transverse  direction,  and  is  the  parallel  component  along  the  axis  of 
undulator.  We  see  that  the  electron  energy,  is  a  constant  of  motion  in 

the  magnetic  field  above.  Substituting  (1)  into  (2)  gives 

B  [-sin  {k  ),cos  {k  ^  ),0]  .  (3) 

at  ^  me 


The  velocity  can  be  solved  exactly  from  (3),  but  assuming  perfect  injection, 
y  :s>  1,  and  (3^  =  1,  the  velocity  components  can  be  written  as 


F(n  = 


K  K 

— cos  ik  (fit ), - sin  {k  (fit  ),P^ 

Y  Y 


where 


1- 


1+K- 


f 


1/2 


(4a) 

(4b) 


and  K  =  eBXQl2jimc’^  is  the  dimensionless  undulator  parameter  or  vector 
potential. 

The  electron  position  can  be  found  by  an  additional  integration  of  (4a) 
and  is 


nt)  = 


KXq  KXq 

-sin  (coo^  ),-7^ — cos  (coof  ),P2 ct 


2ky 


2jrY 


where  coq  =  P^/^o^•  ^  typical  FEL,  /C  =  1,  Xq  =  5  cm,  andy  =  100,  so  that 

the  electrons  travel  at  speed  p^c  =c  along  the  z-axis  for  several  meters 
while  executing  small  transverse  oscillations  with  amplitude  KXol2ny=  0.1 
mm,  which  is  smaller  than  the  typical  beam  radius.  The  transverse 
deflection  of  the  highly  relativistic  electron  beam  in  the  alternating,  periodic 
field  of  the  undulator  leads  to  spontaneous  emission.  Even  if  the  deflection 
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is  too  small  to  be  visualized  on  the  scale  of  the  undulator,  it  is  still 
responsible  for  crucial  spontaneous  emission,  and  FEL  gain. 

The  spontaneous  emission  mechanism  is  essentially  classical,  and  can 
be  described  by  the  Larmor  formula  giving  the  total  power  emitted  from  a 
single  electron  as  Pg  =  2c^P^/3c  ~  2'^B^cr^l3,  where  =  e'^lmc^  is  the 
classical  electron  radius  (in  cgs  units)  [12].  This  radiation  is  emitted  into  a 
small  cone  about  the  forw-ard  direction  with  angular  width  =  7^^.  This 
narrow -spectrum  radiation  is  trapped  in  the  optical  resonator  and  used  to 
supply  the  starting  photons  for  oscillator  operation.  The  interaction 
between  these  photons  with  the  electron  beam  in  the  undulator  is  the 
essential  concept  of  FEL  physics.  In  the  amplifier  configuration,  the 
photons  are  supplied  by  the  seed  laser. 

The  fundamental  feature  of  the  FEL  mechanism  is  that  the  beam 
electrons  undergo  axial  bunching  in  the  combined  undulator  and  radiation 
fields  to  exchange  energ>’  with  an  optical  field.  Assuming  a  circularly 
polarized  plane  wave,  the  optical  fields  are  taken  to  be  [12] 

=  E(cosv,-sinv,0)  ,  and  =  E  (sin  \j/,cos  V|/,0)  ,  (5) 

where  E  is  the  optical  field  strength,  =  kz  -cof  +(^  with  carrier  frequency 
O)  =  kc ,  k  =  2n/X  is  the  optical  w'avenumber  correlated  with  the  wavelength 
X,  and  p  is  the  optical  phase. 

Since  the  electron  travels  at  slightly  less  than  the  speed  of  light,  several 
optical  wavelengths  pass  over  an  electron  as  it  travels  the  full  length  of  the 
undulator.  When  exactly  one  optical  wavelength  passes  over  an  electron  as 
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that  electron  travels  through  one  unduiator  wavelength,  resonance  is 
established.  Efficient  energy  exchange  requires  that  the  electrons 
experience  near  resonant  forces  from  the  optical  and  unduiator  fields.  For 
>  »  1,  this  resonance  is  achieved  when  the  optical  wavelength,  the  beam 
energy,  the  unduiator  strength,  and  the  unduiator  period  approximately 
satisfy  [1] 
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The  expression  for  the  resonant  FEL  wavelength  above  shows  one  of  the 
FEL’s  most  important  attributes,  continuous  tunability.  As  the  electron 
energy  from  the  accelerator  is  varied,  or  as  the  unduiator  strength  K  is 
varied,  the  resonant  wavelength  is  changed. 


B.  OPTICAL  FIELD  EVOLUTION  AND  THE  WAVE  EQUATION 


In  order  to  follow  the  growth  of  the  optical  field  it  is  necessary  to 
develop  a  wave  equation.  Start  with  Maxwell’s  wave  equation  in  the 
Coulomb  gauge  [12], 


where  A  is  the  optical  vector  potential,  and  is  the  transverse  current 
density.  An  appropriate  form  of  the  vector  potential  for  a  circularly 
polarized  optical  wave  is  A*  =  (E/^  )[sinv|r,cosv|/,0].  The  total  current  can  be 
represented  as  the  sum  of  the  single-particle  currents  in  the  electron  beam, 
= -ec ),  where  7^  follows  the  trajectory  of  the  particle. 
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Assuming  that  the  amplitude  and  phase  of  the  optical  field  is  taken  to  vary 
slowly  over  an  optical  period  {E  <*:  oiE ,  ^  c  ox}))  and  over  an  optical 
wavelength  {E'  c  kE ,  a  ki;)  and  keeping  only  terms  with  single 
derivatives,  the  equation  can  be  written  as 

~  -  ,  where  C  =  OiQ+k  )z  -  co^  , 

Y 

p  is  the  electron  beam  particle  density,  the  brackets,  <..>,  indicate  a  phase 
average  over  the  sample  electrons  within  a  wavelength  of  light,  and  ^  is  a 
dimensionless  variable  that  describes  the  electron’s  phase  vdth  respect  to 
the  combined  optical  and  undulator  fields  within  roughly  one  optical 
wavelength,  since  /j  :s>  /jq  2  =  If  we  assume  a  long  electron  pulse 

and  don’t  consider  longitudinal  optical  modes,  we  can  remove  z  -dependence 
from  the  wave  equation,  then  the  wave  equation  can  be  written  as 

a  =  -j<  >  ,  for  0  <  t  <  1  ,  (6) 


A.  +  LA 

dz  C  dt 


where  the  complex  dimensionless  field  envelope  and  current  density  are 


a  = 


AneKLNEe^^ 


J  = 


_  m{eTiKL)\F 


■^mc' 


F  =  r^/w'^  is  the  "filling  factor"  which  is  the  ratio  of  the  electron  beam  area, 


Tir^,  to  the  optical  mode  area  averaged  over  the  undulator  length. 


x-ctiL  is  the  dimensionless  time,  and  the  "  "  denotes  a  time  derivative 
with  respect  to  x.  If  many  electrons  are  randomly  spread  in  longitudinal 
position,  then  the  phase  average  <e"'^>  =  0,  with  no  optical  field  growth. 
In  order  to  understand  the  FEL  interaction,  first  we  must  examine  how 
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electrons  evolve  from  their  initial  positions  within  a  wavelength  of  light. 


C.  ELECTRON  DYNAMICS  AND  THE  PENDULUM  EQUATION 

The  interaction  of  the  electrons  and  the  optical  fields  inside  the 
undulator  is  described  by  the  complete  Lorentz  force  law  and  the  energy 
equation  [12]; 


— (Y^  =  +  (]?x5^)]  ,  (7a) 

at  me 

and  .  (7b) 

at  me 

Inserting  the  full  combined  fields  of  the  undulator  field  (1)  and  optical  field 
(5),  the  transverse  motion  with  perfect  injection  can  be  fovmd  to  be 


K  A 

= - [cos  ik  qZ  ),sin  (k  o2  ),0]  +  — [sin  \)/, 


rns  Ai/Ol 


(8) 


where  A  =eE}J2Kmc^  is  related  to  the  optical  vector  potential  and  is 
analogous  to  the  definition  of  K. 

From  (7b)  and  (8),  the  change  of  electron’s  energy  is  found  to  be  [1] 

_  eEK^  ^ 

yme 

The  exchange  of  energy  of  the  electrons  is  easily  xmderstood  in  terms  of  the 
electron’s  phase,  For  example,  electrons  with  phases  such  that, 
-k/2  <  (^+4>)  <  k/2,  will  get  energy  from  the  optical  field  causing  them  to 
accelerate.  Other  electrons  with  phases,  k/2  <  <  3k/2,  will  lose  energy 

to  the  optical  wave  causing  it  to  grow.  For  a  random  distribution  of  initial 
phases,  about  half  the  electrons  gain  energy  and  about  half  lose.  But,  under 
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the  influence  of  the  combined  fields,  the  electrons  will  begin  to  "bunch "  at 
desirable  phases  and  cause  gain  during  a  pass  down  the  undulator. 

To  further  clarify  the  dynamics  of  the  electrons  in  the  undulator,  it  is 
convenient  to  rewrite  (9)  in  terms  of  ^  alone.  Using  the  relation 
and  the  conditions  1,  P^  =  1,  and  P^  it  can  be 

shown  that  the  individual  electrons  follow  motion  described  by  a  simple 
pendulum  equation  [1] 

00  o 

c  =  V  =  lalcos(C+<l))  ,  (10) 

O 

where  v  =  ^  =  L  [(/j  o^Pa  ]  is  the  dimensionless  phase  velocity. 

The  wave  equation  (6)  and  the  pendulum  equation  (10)  are  valid  in  both 
strong  and  weak  optical  fields  with  high  or  low  gain.  They  are  not  accurate 
for  the  high-efficiency  FEL  where  there  are  significant  changes  in  the 
electron  energy  ymc  The  electron  energy  is  treated  as  a  constant  along  the 
undulator  in  the  definition  of  the  dimensionless  field  strength  la  I,  and  is  not 
followed  self-consistently.  The  optical  field  strength  la  I  explains  why 
stronger  optical  fields  produce  rapid  electron  evolution  and,  how  more 
relativistic  electrons,  with  a  larger  y,  require  a  stronger  optical  field  E  to 
accomplish  the  same  bunching.  When  an  electron  loses  energy,  the  value  of 
la  I  increases  so  that  the  coupling  between  ^  and  v  increases.  When  more 
than  10%  of  the  electron  beam  energy  is  converted  to  radiation,  the  FEL 
dynamics  can  be  altered  so  that  the  more  complete  wave  and  pendulum 
equations  may  be  used. 
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III.  BEAM  QUALITY  AND  GAIN  EVOLUTION 


A  THE  FEL  INTEGRAL  EQUATION 

To  describe  the  FEL  classically,  the  electron  dynamics  are  governed  by 
the  optical  wave  (6)  and  electron  pendulum  (10)  equations.  In  the  case  of 
weak  fields,  la  1  «  Ji,  FEL  dynamics  can  be  expressed  in  the  form  of  an 
integral  equation  by  solving  the  electron  pendulum  and  optical  wave 
equations  self-consistently  [1,  5-7], 

a(x)  =  f^dx'  X'  F(t')  e a (x-xO  ,  (11) 

where  F(x')  =  jdqfiq)  e"**?*^  is  the  characteristic  function  of  the  distribution 
f(q),  f(q)  is  the  distribution  of  initial  electron  phase  velocities  Vj  =VQ+q 
about  Vq,  and  jf{q)dq  -  1.  Note  that  F(x')  and  fiq)  only  depend  on  the 
electron  beam’s  initial  conditions  with  subsequent  dynamics  described 
exactly.  The  FEL  integral  equation  (11)  is  generally  solved  numerically,  and 
has  a  distinct  advantage  over  simulations  which  must  follow  many  tens  of 
thousands  of  sampled  electrons. 

When  the  current  is  low,  j  <  n,  the  field  a  (x-x')  can  be  approximated  by 
the  initial  field  a(0)  =  ao  in  the  integrand  of  (11).  In  the  case  of  perfect 
beam  quality,  f  iq)-  6(0)  and  F(x')  =  1,  the  field  at  the  end  of  the  undulator 
a  (1)  is  easily  evaluated  to  find  the  spectrum  for  low-gain, 
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G  (vq)  -j  [2  -  2cos  (vq)  -  Vosin  (vq)]  /  (12) 

where  gain  is  defined  as  G  =  (a^(l)  -  a^)  I  a  q  .  The  antisymmetric  shape  of 
the  low-gain  spectrum  has  width  Avq  =  2k  and  a  peak  value  G  =  0.13J  at 
Vq  =  2.6.  These  features  are  among  the  first  FEL  characteristics  found 
experimentally  [1]. 

When  the  current  is  high,  J  :>  n,  the  optical  field  grows  exponentially 
and  at  resonance,  Vq  =  0,  where  the  growth  rate  is  maximum,  the  gain  is 
G(t)  =  exp  ((j'72)^''^\'3x)  /  9.  Near  resonance,  Vq  =  0,  the  final  gain  at  t  =  1  is 

G  (vq)  =  exp(V30  /2)^'^)  (  1  -  /  9  •  (13) 

The  gain  spectrum  for  high  current  is  symmetric  about  resonance,  and 
drops  off  at  a  characteristic  value  of  the  phase  velocity 
Avq  =  =  4.22j^'^.  For  large  j,  say  j  =  30,  the  peak  gain, 

G  =  exp(V3072)^'3)/9  =  8  is  much  larger  over  a  broader  spectrum,  Avq  =  7.5, 
than  in  the  low-gain  case. 

B.  ELECTRON  DISTRIBUTION  FUNCTIONS 

Many  experiments  do  not  have  the  perfect  beam  quality  used  in  the 
anal5rtic  results  (12)  and  (13).  Depending  on  specific  experimental 
conditions,  the  actual  shape  of  the  electron  beam  distribution  function  may 
vary  from  one  experiment  to  another.  An  initial  electron  phase  velocity 
Vj  =  L[(k  -t- ^o)Pz((^)  ~  ^  1  depends  on  the  initial  electron  z  velocity  cP^(O), 
and  is  therefore  a  function  of  both  the  energy  spread  and  angular  spread  in 
the  beam.  Consider  a  relativnstic  electron  beam  with  average  energy  ymc^, 
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and  spread  described  by  width  Oq  =  4tiN Ay/y.  The  resulting  distribution  to 
be  used  in  the  integral  equation  (11)  is  the  Gaussian  distribution  [5-7], 

-9*/2a3 

f(q)  =  -^= -  for  all  q  .  (14) 

Alternatively,  an  angular  spread  caused  by  emittance  in  the  beam  can 
result  in  the  exponential  distribution  for  f  iq).  A  relativistic  electron 
entering  the  undulator  with  velocity  at  a  small  angle  0  has  a  reduced  z 
velocity,  =  cPqCOsG  =  cPo(l~^^''2),  and  the  corresponding  change  in  the 
phase  velocity,  v(0)  =  Vq  -  2kN^Q‘^/(.1+K^).  A  Gaussian  spread  in  angles  0  of 
width  A0  results  in  the  exponential  distribution  [5-7], 

f iq)  = -  for  q  <0  and  f iq)  =  0  for  q>0  ,  (15) 

where  Cq  =  4TCA''7^A0^/(l-i-/!r^).  If  a  very  broad  angular  distribution  is  cutoff 
by  an  emittance  filter,  the  resulting  phase  velocity  spread  is  a  triangular¬ 
shaped  distribution  with  positive  slope  [7], 

/■  (q )  = -^(a+ +  q )  for  -c^.<q<0  and  f{q)  =  0  otherwise  ,  (16) 

where  o+  is  the  width  of  the  triangle. 

In  a  realistic  FEL,  the  correct  shape  of  fiq)  is  generally  not  knowm 
precisely,  and  is  determined  by  several  complicating  factors  involving  the 
accelerator  and  electron  beam  transport  properties.  The  L  =  5m  Paladin 
experiments  are  taken  to  be  described  by  the  Gaussian  distribution  (14), 
whereas  the  L  =  15m  Paladin  experiments  made  use  of  an  emittance  filter 
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and  transport  system  that  resulted  in  a  distribution  more  like  (16)  or 
possibly  (15). 

C.  GAIN  SPECTRUM  BROADENING  IN  THE  5  m  PALADIN 

As  the  LLNL  Paladin  experiment  developed,  measurements  were  made 
with  undulators  of  first  5  m  and  then  later  15  m  length  [2-4].  In  the  5  m 
experiment,  the  electron  beam  current  was  about  600  A  at  about  ymc^  =  41 
MeV  energy  and  had  a  radius  of  -  5  mm.  The  L  =  NXq  =  5  m  undulator 
contained  N  ~  62  periods  of  wavelength  Xq  =  8  cm,  and  undulator  parameter 
K  =  0.9.  The  radiation  wavelength  was  X  =  10.6  pm  in  a  radiation  mode 
with  waist  Wq=  4  mm  that  was  slightly  smaller  than  the  electron  beam 
radius  ~  5  mm  at  the  beginning  of  the  undulator  so  that  it  reduced  the 
effective  current  density.  Natural  diffraction  allows  the  optical  mode  to 
expand,  and  we  determine  the  filling  factor  F  in  J  by  averaging  along  the 
undulator.  As  a  result,  we  describe  the  5m  Paladin  experiment  by 
dimensionless  current  j  =  30.  The  observed  electron  beam  brightness 
establishes  that  the  spread  in  initial  phase  velocities  is  about  Avq  =  12,  but 
the  shape  of  the  distribution  is  not  known.  We  take  the  shape  to  be  the 
Gaussian  distribution  as  described  in  (14)  with  Cq  =  12. 

In  figure  1,  the  weak-field  gain  spectrum  (shown  as  dots)  was  measured 
by  varying  the  undulator  field  strength  B  over  a  small  range  to  change  the 
resonance  condition  Vq  [5-7].  The  observed  gain  spectrum  width  is  Avq  =  40, 
with  peak  gain  of  G  =  12%.  The  solid  curve  uses  the  integral  equation  to 
evaluate  gain  with  j  =  30  and  Cq  =  12.  The  theoretical  gain  spectrum  was 
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observed  to  be  antisymmetric  in  shape  with  a  peak  gain  of  G  =  0.15.  The 
observed  and  theoretical  gain  spectra  are  dramatically  distorted  by  the 
effects  of  poor  beam  quality.  Figure  2  shows,  for  comparison,  the  gain 
spectrum  determined  by  the  integral  equation  with  j  =  30  and  =  0 
describing  perfect  beam  quality.  Poor  beam  quality  is  responsible  for 
decreasing  the  FEL  peak  gain  from  G  =  9  in  figure  2  to  G  =  15%  in  figure  1, 
and  distorting  the  gain  spectrum  from  an  approximately  symmetric  shape  of 
width  Avq  =  6  to  an  antisymmetric  shape  of  width  Avq  =  70.  Note  that  the 
observed  and  theoretical  antisymmetric  spectra  in  figure  1  are  both  much 
broader  than  the  antisymmetric  low-gain  spectral  shape  in  (12)  with  width 
Avq  =  2k. 

The  shape,  overall  width,  and  peak  value  of  the  5m  Paladin  gain 
spectrum  is  dramatically  altered  from  the  idealized  gain  spectrum  (either 
(12)  or  in  figure  2),  and  in  this  case,  is  determined  by  the  initial  distribution 
of  electron  velocities  (14). 
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Figure  1.  Gain  evolution  of  the  5  m  Paladin  with  poor  beam  quality 
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D.  A  PLATEAU  IN  THE  GAIN  EVOLUTION  OF  THE  15  m  PALADIN 


We  also  identify  characteristics  of  the  Paladin  15m  experiments  that 
are  imique  to  the  initial  electron  distribution  function  shape.  An  emittance 
filter  was  used  to  improve  the  beam  quality  entering  the  undulator,  and  a 
detailed  study  of  the  filter  and  transport  system  concluded  [13]  that  the 
distribution  function  shape  should  look  something  like  the  exponential  or 
triangular  distributions,  (15)  or  (16). 

In  the  15  m  experiment,  the  electron  beam  current  was  about  500  A  at 
ymc^  =  45  MeV  energy.  The  longer  undulator  was  L  =  15  m  long, 

containing  N  ~  187  periods  with  K  =  1.0.  The  high  gain  observed  in  the  15 
m  experiment  causes  the  optical  mode  and  electron  beam  to  be  well  matched 
over  the  whole  undulator  length,  so  that  we  take  F  =  1.  Other  parameters 
are  the  same  as  in  the  5  m  experiment.  Because  of  the  longer  undulator 
length,  the  15  m  experiment  is  described  by  large  j  =  1400,  and  the 
triangular  distribution  (16)  with  =  13.  For  such  a  large  current,  the  ideal 
gain,  with  perfect  beam  quality,  a+  =  0,  would  be  G  =  5x10^  peaked  at  Vq  =  0 
and  would  evolve  exponentially  in  t.  Figure  3  shows  the  theoretical  gain 
evolution,  G(t),  found  by  solving  the  integral  equation  with  the  triangular¬ 
shaped  distribution  (16)  of  width  =  13,  j  =  1400,  and  Vq  =  0.  The  gain 
evolution,  G(t),  represents  the  optical  power  growth  along  the  undulator 
length  as  t  goes  from  0  to  1,  and  it  shows  a  dramatic  plateau  about  half  way 
along  the  undulator  from  x  =  0.5  to  0.7.  The  choice  of  an  exponential 
distribution  (15)  shows  significantly  smaller  gain  than  the  triangular-shaped 
distribution  and  with  a  more  subtle  plateau. 
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Figure  4  shows  the  optical  power,  P,  measured  along  the  15  m 
undulator  length  [2],  corresponding  to  x  =  0  ->  1.  The  changes  in  power 
indicate  the  gain  evolution,  and  the  figure  shows  a  the  prominent  plateau  in 
the  middle  of  interaction  length.  It  appears  that  the  shape  of  the 
distribution  function  in  the  Paladin  15  m  experiment  has  resulted  in  a 
plateau  in  the  FEL  interaction  half  way  through  the  undulator,  and  the 
experimental  plateau  in  the  gain  evolution  is  verified  by  the  integral 
equation  with  the  positive  sloped  triangular-shaped  distribution. 

In  conclusion,  the  results  of  the  Paladin  5  m  and  15  m  experiments, 
have  provided  a  unique  opportunity  to  observe  how  the  FEL  gain  spectrum 
and  gain  evolution  can  be  determined  by  the  electron  beam  velocity 
distribution  instead  of  the  natural  high-gain  mechanism.  The 
measurements  provide  the  first  comparison  between  theory  and  experiment 
in  a  short-wavelength  device  where  the  electron  beam  distribution  plays  a 
dominant  role  in  the  FEL  mechanism.  The  gain  spectrum  shape,  width, 
and  evolution  can  be  used  to  determine  unique  information  about  beam 
quality. 
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Figure  3.  Theoretical  gain  along  the  15  m  Paladin  undulator 
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IV.  A  SIMPLE  MODEL  OF  THE  ELF  I  EL 


A  DESCRIPllON  OF  ELF  EXPERIMENT 

The  LLNL  Electron  Laser  Facility  (ELF)  FEL  is  an  amplifier  with  very 
high  gain  operating  in  the  microwave  regime  using  the  Experimental  Test 
Accelerator  (ETA).  The  experiment  examined  high-power  saturation,  and 
the  effect  of  tapering  the  undulator  field.  ETA  uses  a  6  kA  electron  beam 
with  a  normalized  emittance  of  I.Stc  cm-rad  [8].  An  emittance  filter  reduces 
the  beam  current  to  approximately  450  A  over  a  pulse  length  of  15  ns  wdth 
normalized  edge  emittance  of  0.4771  cm-rad.  The  beam  energy  is  3.3  MeV 
with  Y  =  7.5.  The  electron  beam  particle  density  is  p  =  2  x  10^^  cm"^,  and 
the  beam  radius  is  =4  mm. 

The  L  =  3  m  long  ELF  undulator  is  linearly  polarized,  with  wavelength 
Xq  =  9.8  cm  and  undulator  constant  K  ~  2.8.  The  peak  field  strength  is 
B  ~  4300  G  over  N  =  30  undulator  periods.  The  Rayleigh  length  is  less 
than  1%  of  the  undulator  length  so  that  ELF  must  use  a  3  cm  x  10  cm 
rectangular  wavegmde  to  confine  the  radiation  wavefront  around  the 
electron  beam  for  continued  amplification.  Without  a  waveguide,  natural 
diffraction  would  spread  the  optical  wave  cross-section  away  from  the 
amplifying  electron  beam  and  reduce  coupling.  The  waveguide  cutoff 
frequency  is  5  GHz  in  the  TEqi  mode,  so  that  the  FEL  at  35  GHz  operates 
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far  above  cutoff. 


B.  SIMPLE  FEL  EQUATIONS  FOR  ELF  EXPERIMENT 

The  FEL  equations  described  in  previous  chapters  are  the  self- 
consistent  pendulum  and  wave  equations.  They  provide  a  simple  model  of 
the  FEL  that  is  valid  in  strong  and  weak  optical  fields  for  high  and  low 
beam  current.  The  assumptions  required  to  derive  the  simple  form  of  the 
equations  are  y  ^  1  and  N  1.  Since  the  most  important  aspects  of  the 
interaction  take  place  near  resonance,  the  changes  in  the  electron  energy, 
Ay/y,  are  small  when  iV  »  1,  and  the  FEL  efficiency  is  low.  For  high 
efficiency,  in  the  tapered  undulator  or  in  a  short  undulator  with  few  periods, 
the  changes  in  the  electron  energy  are  expected  to  be  large.  High  efficiency 
FELs  the  undulator  can  be  tapered  to  maintain  resonance  w’ith  the  trapped 
electrons  that  lose  energy  to  the  optical  wave.  The  untrapped  electrons  will 
drift  far  off-resonance.  Tapering  and  high-efficiency  can  be  added  to  the 
basic  FEL  electron  equation  of  motion  (10),  and  it  becomes  [14,  15] 

oo  o 

;  =V  =  e(lal -aJ5  +  (l-(3v/4aN))lal  cosCC  +  b)  ,  (17) 

where  Q(x )  is  the  step  function,  Q(x )  =  0  for  x  <  0  and  0(a:  )  =  1  for  x  >  0,  so 
that  the  tapering  starts  at  the  onset  of  saturation,  Og  =  2{j  /  2)^^  is  the 
saturation  field,  5  =  4kN Ay^  /  y^  is  the  taper  rate  (5  =  0  for  the  untapered 
undulator),  and  y^  is  the  resonant  Lorentz  factor.  The  modified  pendulum 
equation  (17)  relaxes  the  assumption  of  N  3>  1,  but  still  assumes  y:»  1  and 
allow  tapering  of  the  undulator. 


23 


In  order  to  have  a  self-consistent  theory,  the  dimensionless  single-mode 
wave  equation  is  also  altered  to  account  for  high  efficiency  in  terms  of  the 
electron  phase  and  phase  velocity  [14,  15], 

a  = -7  <  (1  -  (v  /  47dV))  >  .  (18) 

In  these  formulations,  electrons  near  resonance,  v  =  0,  are  described 
accurately,  while  imtrapped  electrons  far  from  resonance,  v  — >  4nN ,  are 
handled  less  accurately.  However,  the  untrapped  electrons  are  expendable 
since  they  become  randomly  spread  in  phase  as  well  as  uncoupled  from  the 
bunching  interaction.  For  electrons  near  resonance,  the  bunching  rate  is 
determined  by  the  field  strength  la  I.  When  lo  I  <t;  k,  the  optical  field  is 
considered  weak,  and  bunching  is  small.  When  la  I  »  Ji,  i.e.,  the  optical  field 
is  strong,  the  electrons  bunch  rapidly,  and  become  trapped  in  closed  phase- 
space  orbits.  The  current  j  determines  the  couph’ng  between  the  electron 
beam  and  the  optical  wave.  Without  current,  j  =  0,  or  no  bunching  <..>  =  0, 
there  is  no  change  in  the  initial  optical  field,  a(0)  =  ao.  The  modified 
pendulum  equation  (17)  and  wave  equation  (18)  are  valid  for  weak  fields, 
strong  fields,  low  gain,  high  gain,  tapered,  untapered,  low  efficiency,  and 
high  efficiency,  but  assume  y  ^  1- 

In  a  linearly-polarized  imdulator,  each  electron  experiences  fast, 
periodic  oscillations  that  are  comparable  to  the  radiation  wavelength  and 
modify  the  interaction  strength.  For  this  coupling  [1],  the  dimensionless 
current  density  j  should  include  the  factor  (e7o(^)  -  e/i(v)^  =  0.54  where 
^  /  2(1  -f-  K^)  =  0.44  for  ELF.  For  early  ELF  experiments  with  450  A 
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beam  current  [9],  the  filling  factor  is  F  =  0.026,  and  the  corrected  current  is 
j  =  3015. 

A  relativistic  electron  in  a  rectangular  waveguide  far  above  cutoff  has 
its  phase  velocity  shifted  by  [1] 


Av 


TtATXo^d  4-^2) 

8y^ 


(19) 


where  X  is  the  long  waveguide  dimension,  Y  is  the  short  waveguide 
dimension,  and  p  and  q  are  integers.  For  ELF  in  the  mode,  the 

resonant  phase  velocity  is  shifted  by  Avqi  =  -  6;r.  A  change  in  the  imdulator 
field,  6B,  about  the  resonant  field  strength  Bq  also  causes  a  change  in 
phase  velocity  by 


4kNK^  M 
(1  +  ^2)  B 


Resonance  in  the  ELF  FEL  is  tuned  using  changes  in  the  magnetic  field 
dB  /  Bq  [8].  Combining  (19)  and  (20)  shows  that  the  waveguide  increases 
the  resonant  undulator  field  strength  by  ABoj  =  240  G. 

The  distribution  of  initial  phase  velocities,  /"(v),  is  important  to  the  FEL 
growth  rate  and  in  ELF  is  determined  by  the  electron  beam  emittance  filter. 
The  emittance  filter  results  in  a  triangular-shaped  distribution  function 
with  a  positive  slope  given  by 

/■+(v)  =  (v  -  Vo  +  2a+)  /  20+  ,  for  Vq  -  <  v  <  Vq  , 

where  Vq  is  the  peak  of  the  beam  distribution.  The  full-vvddth  at  the  base  of 
the  triangle  is  equivalent  to  a  6.4%  energy  spread  so  that  2a^  = 
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4kN  X  0.064,  and  o+  =  12  [10].  The  initial  input  power  in  ELF  is  30  kW  in 
the  TEqi  waveguide  mode,  and  corresponds  to  a  dimensionless  microwave 
field  magnitude  of  Uq  =  12. 

C.  RESULTS 

Figure  5  shows  the  simulated  gain  spectrum,  G(vo),  and  figure  6  shows 
the  experimental  result,  G(Avg),  at  the  end  of  the  undulator.  Recall  that 
the  equation  (20)  shows  the  relationship  between  the  change  in  Vq  used  in 
the  simulation  and  the  change  in  the  magnetic  field  used  in  the  experiment. 
The  data  is  taken  with  I  =  450  A  beam  current  so  that  j  =  3015,  and  = 
12  [9].  The  simulation  gives  peak  gain  Gn,ax  ~  2833  at  Avq  =  28  away  from 
resonance  with  a  gain  bandwidth  of  6vo  =  30.  The  observed  peak  gain  is 
four  times  higher  than  the  calculated  result,  and  is  shifted  from  resonance 
more  than  in  the  simulation.  However,  the  important  qualitative  features  of 
the  experiment  are  shown  by  the  simple  single-mode  theory.  The 
asymmetrical  gain  spectrum  shape  is  in  good  agreement.  Mode  distortion 
by  excitation  of  higher-order  waveguide  modes  probably  accounts  for  the 
some  of  higher  measured  gain,  and  the  larger  shift  away  from  resonance  in 
the  experiment.  Some  of  the  disagreement  may  also  be  due  to  experimental 
error,  since  the  powerful  microwave  pulse  must  be  attenuated  by  many 
orders  of  magnitude  before  a  measurement  can  be  made. 
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The  extraction  of  higher-order  modes  has  the  effect  of  confining  the 
microwave  power  around  the  electron  beam,  increasing  the  effective  filling 
factor,  and  increasing  coupling  in  the  later  parts  of  the  undulator.  In  order 
to  account  for  mode  distortion,  we  increase  the  filling  factor  from  0.026  to 
0.032.  Figure  7  shows  the  simple  simulation  results  for  power  evolution, 
P(x),  and  figure  8  shows  the  experiment  results  for  power  evolution,  P(2). 
The  later  ELF  experiments  increase  the  current  to  850  A  [10,  11]  so  that 
the  dimensionless  current  density  with  the  new  filling  factor  is  now  j  = 
7000.  The  dimensionless  microwave  field  is  Oq  =  15,  and  =  12.  The 
onset  of  saturation  occurs  at  x  =  0.5,  or  about  1.5  m  along  the  undulator, 
and  the  efficiency  is  t]  =  5%  in  the  untapered  case.  The  shape  of  P(x)  in  the 
untapered  case  in  figure  7  and  figure  8  shows  excellent  agreement  between 
the  simple  theory  and  experiment. 

In  order  to  extend  the  natural  saturation  limit,  the  undulator  properties 
can  be  altered  along  z  to  maintain  resonance  and  interaction  strength.  The 
phase  acceleration,  5,  is  given  by  h~-AnNK‘^/SB  when  the 

undulator  field  strength  is  decreased  by  AP  [16].  Figure  7  shows  the  result 
of  tapering  with  5  =  90k  giving  an  efficiency  of  ri  =  31%.  The  experimental 
results  in  figure  8  shows  that  the  efficiency  is  increased  from  6%  to  35%, 
because  of  tapering  [11].  The  efficiency  is  increased  by  about  a  factor  of  six 
in  both  simulation  and  the  experiment.  In  both  cases,  tapered  and 
untapered,  we  see  excellent  agreement  between  the  experiment  and  the 
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simplified  simulation  . 


In  each  comparison  above,  the  excitation  of  higher-order  waveguide 
modes  confines  the  microwave  field  tightly  around  the  electron  so  that  the 
filling  factor  F  and  the  coupling  are  effectively  increased.  The  ELF 
experiments  determined  that  about  30%  of  the  microwave  power  is  in 
higher-order  modes  [10],  consistent  with  this  model. 
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Figure  8.  Gain  evolution  of  ELF  experiment 
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